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J^BSTRACT 

TMs routine generates a (pseudo) random sequence of integers 
in the range ^ N 4 2**16-1, 
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1- RSQUI RE GENTS 
!•! Memory 

IK or larger alterable memory 

1.2 Equ.lpment 

NOVA central processor 

1 . 3 Externa 1 Svhrgntine s 
None 

1.4 Other 
None 

2. OPEimTING PROCEDUR g 

2 ^ 1 C alli n g secmence 

JSR . FJ^ND 

address of old value 

return 

2.2 Input Format 

The address of the previous random value (or initially 
a starting value) must be provided in the word after 
the JSR .RAND. 

2.3 Output For mat 

The new 16-bit random result will be returned in IhC0 
and will also replace the previous value in memory. 

2.4 Error Returns 
None 

2.5 state of l^ctive Registers upon Exit 

ACl and AC2 are unchanged. hC0, AC3 and Carry are 
destroyed. 
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2.6 Cautions to User 

If a K-bit nun&er (1 ^ K <16) is needed^ use the most 
significant K bits (the least significant K bits are 
not as random) • For example, to obtain random N mod 2, 
use the sign bit of the result, 

3. DISCUSSION 

3.1 Algorithms 

This routine generates a linear, cozigruential sequence 
of the form 

X(N + 1) = (X(N) * A + C) mod 2**16 

where N = (2^*11 + 2*^*2 + 1) 
and C = 33031 (octal) 

3.2 Limitations and Accuracy 

The resulting sequence has the maximum period of 
2'**16. 

3.3 Size and Ti! T>ing 

The routine is 36 (octal) words in length. 
Execution time is 244.7 /t seconds. 

3.4 References 

Donald E. Knuth, The ^xt of Computer Programming , 
Volume 2, Addison-Wesley Publishing Company (1969), 
Chapter 3 . 

3.5 Flow Diagrams 
None 

4. EXAMPLES AND ^PPLICATIOMS 

Although this routine produces numbers in the range 

4 N $ 2**16 - 1 
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fractional values in the range 

^4 F < 1 
can be obtained by tlie simple computation 

F «- N/2**16 

Random nuntoers are useful in a vBiletY of applications - 
simulation of natural phenomena, random sampling, testing 
the effectiveness of computer algorithms, and probability 
theory to name a fei-f. 

The source of .RAND is pro¥ided with the NOV^ software. 
If a user routine requires this program, the tape should 
be edited into the user source. 

5. PROGRAM LISTING 

A listing of .RAND follows. No origin is given in the 
source, enabling the user to edit this tape anywhere within 
his routines. 
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WORD 



OUTPUTS 



A-'^s ;.■•» S' Cf^pf-^ ^E°LAC^v^ 3L=3 VALUE 



CALLING SE0UENCE5 

JSR •RAND 
ADDRESS OF OLD VALUE 
RETURN 



METHODS 



CAUTION! 






'-l"'- '. ">''.'^EI5' A Ll'\£^"' V.'J-i^P.'jl'l'SZl :'\l. 



V f .\ .-. . ^ := < X C : .■ > '< 4 '- C 5 ?'' D 2 -'' "- ' ': 
— A x-BiT '<^';3C'" .^='/"^^.:- CAS C-^POSFJ} 

V <'-' " H^ *io ST i^ : (-M 1 5 ! G ^.^: T 

;..-—- i-f'-z l::^:^' si ;-<:": ?^?^t k-bits 
^ ... ._. ^, ,-^^ „ ^^^ ^^. ^. ,j^ . -^ „ .^, ^,^., J ^ 

r-;.. ;: f i;f^7r„f , 7''; j'"''A"A- ■^%-x?:'G>i v mod ?# 



ACl^ AC2 

AC0, AC3.f CARRY 



f^9?'0 IS???::?? 



RANDt 



ST A 


3r- r J -I 2 


* 


SAVE RETURN 


J ";•' 


*..Te'3 


I 


BUMP PAST ADDRESS 


•t: r . 


^ ? .- i;3i ■ 


» 


*SAVE AC I 


- T ?> 


2^- .:'■••■ 2 


1 


♦SAVE AC2 


. "■•..-. 


O ' "• . 


1 


(5ET ADDRESS OF OLD 


* T A 


"/■ ~, ^. '• _^ ^ 


1 


OLD VALUE TO AC0 


• --p 


^ L-Dr>C 


1 


M*A 




'> ..,,■;• '■'"•'^ 


I 


GET INCREMEMT* C 


^ !-,-■. 


" ■ '.f 


1 


Cf^*A+C> MOD 2** I 6 


«^T -^ 


2.»c?.«2 


S 


STORE IT 


lDA 


l*^L--i'l 


1 


♦RESTORE ACl 


uD^ 


2.T »UC15?2 


1 


* RE STORE AC2 




?. , • 5 r. :: -• 


1 


RETURN 



CONSTANT 



VALUE 



•- '-: „ 1 O I- 4 -i C ^'f. '■ 
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i»i31 i0«i»S .UDC2S i 

iii32 tsiiii •uoi3s i 

iii33 iiiiii .uoiii i 

i§i34 §33031 .lIDgis 33i31 

iii35 0§§ili .UD21I li 



I *SAIZ A0£ 

I SAVE RETUsM 

I I IE X AT J CM COUNT 

I : T L HAT 1 '^J CO UN T* K 



